function f=groll(V)

global Vw cl_g cdg T_W mibrake W_S fi rho Nn_W mi g

f= (V+Vw)./(g*((T_W-mibrake)-(((cdg-(mibrake*cl_g))/W_S)*((1/2)*rho*(V.^2)))+...
    (Nn_W*(mibrake-mi))-fi));